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g ; Abstract 

Matrix representations of the Maxwell equations are well-known. However, 
all these representations lack an exactness or/and are given in terms of a 
pair of matrix equations. We present a matrix representation of the Maxwell 
equation in presence of sources in a medium with varying permittivity and 
permeability. It is shown that such a representation necessarily requires 8x8 
matrices and an explicit representation for them is presented. 
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1 Introduction 

O 

""^ ■ Matrix representations of the Maxwell equations are very well-known 

However, all these representations lack an exactness or /and are given in 
terms of a pair of matrix equations. Some of these representations are in free 
space. Such a representation is an approximation in a medium with space- 
On 1 and time-dependent permittivity e(r, t) and permeability /i(r, t) respectively. 

Even this approximation is often expressed through a pair of equations using 
3x3 matrices: one for the curl and one for the divergence which occur in 
the Maxwell equations. This practice of writing the divergence condition 
separately is completely avoidable by using 4x4 matrices [0 for Maxwell 
equations in free-space. A single equation using 4x4 matrices is necessary 
and sufficient when e(r,t) and are treated as 'local' constants Jl|, [3J. 

A treatment taking into account the variations of e(r,t) and n(r,t) has 
been presented in [Q. This treatment uses the Riemann-Silberstein vectors, 



F (r,t) to reexpress the Maxwell equations as four equations: two equations 
are for the curl and two are for the divergences and there is mixing in F + (r, t) 
and F~ (r, t). This mixing is very neatly expressed through the two derived 
functions of e(r,t) and fi(r,t). These four equations are then expressed as a 
pair of matrix equations using 6x6 matrices: again one for the curl and one 
for the divergence. Even though this treatment is exact it involves a pair of 
matrix equations. 

Here, we present a treatment which enables us to express the Maxwell 
equations in a single matrix equation instead of a pair of matrix equations. 
Our approach is a logical continuation of the treatment in [[|. We use the 
linear combination of the components of the Riemann-Silberstein vectors, 
F ± (r, t) and the final matrix representation is a single equation using 8x8 
matrices. This representation contains all the four Maxwell equations in 
presence of sources taking into account the spatial and temporal variations 
of the permittivity e(r,t) and the permeability fi(r,t). 

In Section-I we shall summarize the treatment for a homogeneous medium 
and introduce the required functions and notation. In Section-II we shall 
present the matrix representation in an inhomogeneous medium, in presence 
of sources. 

2 Homogeneous Medium 

We shall start with the Maxwell equations [||, |J in an inhomogeneous medium 
with sources, 

V-D{r,t) = p, 
VxH{r,t)-^D(r,t) = J, 

VxE(r,t) + ^B(r,t)=0, 

V-B(r,t) = 0. (1) 

We assume the media to be linear, that is D = eE, and B = fiH, where 
e is the permittivity of the medium and [i is the permeability of the 
medium. In general e = e(r, t) and fi = fi(r, t). In this section we treat them 
as 'local' constants in the various derivations. The magnitude of the velocity 
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of light in the medium is given by v(r,t) = \v(r,t)\ = 1/ Je(r, t)[j,(r, t). In 
vacuum we have, e = 8.85 x 1(T 12 C 2 / N.m 2 and /i = 47r x 10~ 7 N/A 2 . 

One possible way to obtain the required matrix representation is to use 
the Riemann-Silberstein vector |3j given by 



F~(r,t) = ±(y/^E(r,t)-i-=L=B(r,t)\ . (2) 

For any homogeneous medium it is equivalent to use either F + (r, t) or 
F~ (r,t). The two differ by the sign before 'i' and are not the complex 
conjugate of one another. We have not assumed any form for E(r, t) and 
B(r,t). We will be needing both of them in an inhomogeneous medium, to 
be considered in detail in Section-Ill. 

If for a certain medium e(r, t) and /z(r, t) are constants (or can be treated 
as 'local' constants under certain approximations), then the vectors [r, t) 
satisfy 

i^F±(r,t) = ±vV xF±{r,t)--j={iJ) 
V-F±(r,t) = -L(p). (3) 



Thus, by using the Riemann-Silberstein vector it has been possible to reex- 
press the four Maxwell equations (for a medium with constant e and /i) as 
two equations. The first one contains the the two Maxwell equations with 
curl and the second one contains the two Maxwell with divergences. The 
first of the two equations in (^) can be immediately converted into a 3 x 3 
matrix representation. However, this representation does not contain the di- 
vergence conditions (the first and the fourth Maxwell equations) contained 
in the second equation in (|3|). A further compactification is possible only by 
expressing the Maxwell equations in a 4 x 4 matrix representation. To this 
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end, using the components of the Riemann-Silberstein vector, we define, 



-F+ + iF+ 

Ft 
Ft 
Ft + 



*"(r,*) = 



—F~ - \F: 



F~ - iF- 



(4) 



The vectors for the sources are 



W + 



Then we obtain 



J X -\- ljy 

Jz - Vp 
Jz + Vp 
Jx H~ '^-J y 



w~ = 



Jx iJy 

Jz ~Vp 
Jz + Vp 



(5) 



d_ 

dt 
d_ 

dt 



dt 1 J 



W + 

-v{M* • V}* - - W~ 



(6) 



where '*' denotes complex-conjugation and the triplet, M = (M x , M y , M z ) 
is expressed in terms of 
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(7) 



Alternately, we may use the matrix J = —Q. Both differ by a sign. For 
our purpose it is fine to use either Q or J. However, they have a different 
meaning: J is contravariant and Q is covariant; The matrix Q corresponds 
to the Lagrange brackets of classical mechanics and J corresponds to the 
Poisson brackets. An important relation is f2 = J" 1 . The M-matrices are: 



M, = 



M„ 
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-1 



p. 



(8) 



Each of the four Maxwell equations are easily obtained from the matrix 
representation in This is done by taking the sums and differences of 
row-I with row-IV and row-II with row-Ill respectively. The first three give 
the y, x and z components of the curl and the last one gives the divergence 
conditions present in the evolution equation (H). 

It is to be noted that the matrices M are all non-singular and all are 
hermitian. Moreover, they satisfy the usual algebra of the Dirac matrices, 
including, 

M x (3 = -fJM x , 
M y (3 = -PM V , 



Ml 



MM,, 



~MyM X 



M 2 z =l 



iM z , 

M y M z = —M z M y = iM x , 

M Z M X = -M X M Z = iM y . (9) 

Before proceeding further we note the following: The pair (ty^jM) are 
not unique. Different choices of ^ would give rise to different M , such that 
the triplet M continues to to satisfy the algebra of the Dirac matrices in @. 
We have preferred ^ ± ma the the Riemann-Silberstein vector (Q) in [^J. This 
vector has certain advantages over the other possible choices. The Riemann- 
Silberstein vector is well-known in classical electrodynamics and has certain 
interesting properties and uses 0. 

In deriving the above 4x4 matrix representation of the Maxwell equations 
we have ignored the spatial and temporal derivatives of e(r, t) and /i(r, t) in 
the first two of the Maxwell equations. We have treated e and fi as 'local' 
constants. 



3 Inhomogeneous Medium 

In the previous section we wrote the evolution equations for the Riemann- 
Silberstein vector in (|3p, for a medium, treating e(r,t) and n(r,t) as 'local' 
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constants. From these pairs of equations we wrote the matrix form of the 
Maxwell equations. In this section we shall write the exact equations taking 
into account the spatial and temporal variations of e(r,t) and fi(r,t). It 
is very much possible to write the required evolution equations using e(r, t) 
and fj,(r, t). But we shall follow the procedure in || of using the two derived 
laboratory functions 



Velocity Function : v(r, t) 



Resistance Function : h(r, t) 



1 



y/e(r,t)n(r,t) 



(10) 



The function, v(r,t) has the dimensions of velocity and the function, h(r,t) 
has the dimensions of resistance (measured in Ohms). We can equivalently 
use the Conductance Function, K(r,t) = l/h(r,t) = e(r,t)//i(r,t) (measured 
in Ohms -1 or Mhos!) in place of the resistance function, h(r,t). These 
derived functions enable us to understand the dependence of the variations 
more transparently ||. Moreover the derived functions are the ones which 
are measured experimentally. In terms of these functions, e = 1/y/vh and 

\i = \Jh/v. Using these functions the exact equations satisfied by F* 1 (r, t) 
are 

i^F+(r,t) = v{r,t)(v x F + {r,t)) + ^(Vv{r,t)x F + {r,t) 
v(r,t) 



Vh(r, t) x F (r, t) T=V v ( r > t)h(r, t) J 



2h(r) V ^ ' > ^ ' >) ^2 

i v(r, t) . . , i h(r, t) „ , 
+ — (r, t) + - ' [ F- (r, t) 
2v{r,t) v ' ; 2h(r,t) v ' ; 

x—F- (r, t) = -v(r, t) (v x F (r, t)) - - (Vv(r, t) x F (r, t)] 

~^T) ( Vh{r > t] x F+ (r>t) ) ~ ^\M^m)j 

^^F-ir^+'-^F+^t) 
2v(r,t) v ' 2h(r,t) v ' ' 
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V F (r,t) 



2v(r,t) 
1 



2h(r,t) 



Vv(r,t) ■ F (r,t) 
Vh(r,t) ■ F + (r,t) 



+^yv(r,t)h(r,t)p. 



(11) 



where v = || and h = ||. The evolution equations in (|TT1) are exact (for 
a linear media) and the dependence on the variations of e(r, t) and yu(r, t) 
has been neatly expressed through the two derived functions. The coupling 
between F + (r,t) and F (r,t). is via the gradient and time-derivative of 
only one derived function namely, h(r,t) or equivalently hz(r,t). Either of 
these can be used and both are the directly measured quantities. We further 
note that the dependence of the coupling is logarithmic 



1 



h(r,t) 



Vh(r,t) = V{ln(/i(r,t))} 



1 



h(r,t) 



h{r,t) 



d_ 

Of 



{In (h(r,t))}, (12) 



where 'In' is the natural logarithm. 

The coupling can be best summarized by expressing the equations in (|TT|) 
in a (block) matrix form. For this we introduce the following logarithmic 
function 



C(r,t) 



1 



{lln (v(r,t)) + a x \n(h(r,t))} 



where a x is one the triplet of the Pauli matrices 





' 


1 " 




' -i " 




' 1 


" 




0~ x = 


1 







i 







-1 





(13) 



(14) 



Using the above notation the matrix form of the equations in (|TT| ) is 



ill 



d__d_ 
^di dt 



C 



F + (r,t) 
F (r,t) 



t;(r)(T 2 {lV + V£} x 



F A 
F 



(r,t) 
(r,t) 



7 



V2 



v(r,t)h(r,t)J 



{IV - V£} 



F 
F 



(r,t) 
(r,t) 



+ 



^=y/v(r,t)h(r,t)p. 



V2 



(15) 



where the dot-product and the cross-product are to be understood as 

A u + B v 
Cu + Dv 



A 


B ' 




u 




C 


D 




V 




A 


B 




u 




C 


D 


X 


V 





A x u 
C x u 



B x v 
D x v 



(16) 



It is to be noted that the 6x6 matrices in the evolution equations in (|15|) 
are either hermitian or antihermitian. Any dependence on the variations 
of e(r,t) and fi(r,t) is at best 'weak'. We further note, V (In (v (r, t))) = 
— V (In (n(r, t))) and (In (v (r, £))) = — 4 (In (n(r, t))). In some media, 
the coupling may vanish (V/i(r, t) = and h(r,t) = 0) and in the same 
medium the refractive index, n(r,t) = c/v(r,t) may vary (Vn(r, t) ^ 
or/and n(r,t) ^ 0). It may be further possible to use the approximations 
V (In (h{r, t))) « and | (In (h{r, *))) « 0. 

We shall be using the following matrices to express the exact representa- 
tion 



(17) 



where S are the Dirac spin matrices and ex are the matrices used in the Dirac 
equation. Then, 



" cr " 




a ' 




"10" 


cr 


a = 


a 


1 = 


11 





at 



j 

o 



o 
j 



2h(r,t) 
-v(r,t) 



v(r,t) 



I 





\j30ty 












I 















J 









2u(r,t) 



i/3% 

{M-V + S-w} -i/3 (£ • tu) 

(E* • iw) a„ {AT • V + E* -u} 



(18) 
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where 



u(r,t) = ^-i- Vv(r,t) = {In ^(r,t)} = -~V{hm(r,t)} 

^ (r,t) = 2M^t) V/l(r ' t) = ^ V{ln/l(r ' t)} (19) 

The above representation contains thirteen 8x8 matrices! Ten of these 
are hermitian. The exceptional ones are the ones that contain the three 
components of w(r,t), the logarithmic gradient of the resistance function. 
These three matrices are antihermitian. 



4 Concluding Remarks 

We have been able to express the Maxwell equations in a matrix form in a 
medium with varying permittivity e(r, t) and permeability fi(r, t), in presence 
of sources. We have been able to do so using a single equation instead of a 
pair of matrix equations. We have used 8x8 matrices and have been able 
to separate the dependence of the coupling between the upper components 
and the lower components (^~) through the two laboratory functions. 
Moreover the exact matrix representation has an algebraic structure very 
similar to the Dirac equation. We feel that this representation would be 
more suitable for some of the studies related to the photon wave function. 
This representation is the starting point for the exact formalism of Maxwell 
optics This formalism provides a unified treatment of beam-optics 

and polarization. 
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